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Abstract. We use a new version of our numerical model 
for particle propagation in the Galaxy to study radioactive 
secondaries. For evaluation of the production cross sections 
we use the Los Alamos compilation of all available experi- 
mental cross sections together with calculations using the im- 
proved Cascade-Exciton Model code CEM2k. Using the ra- 
dioactive secondary ratios 26 A1/ 27 A1, 36 C1/C1, 54 Mn/Mn, we 
show how the improved cross-section calculations together 
with the new propagation code allow us to better constrain 
the size of the CR halo. 



1 Introduction 

In recent years, new and accurate data have become available 
in CR astrophysics; more CR experiments are planned for 
launch in several years that will tremendously increase the 
quality and accuracy of CR data making further progress de- 
pendent on detailed models. Data will continue to flow from 
the high resolution detectors on Ulysses, Advanced Compo- 
sition Explorer (ACE) and Voyager space missions. Sev- 
eral high resolution space experiments will be in orbit in 
the nearest 2-3 years, e.g., PAMELA to measure antipro- 
tons, positrons, electrons, and isotopes H through C over the 
energy range of 0. 1 to 200 Ge V, and Alpha Magnetic Spec- 
trometer (AMS) to measure particle and nuclear spectra to 
TeV enegies. 

Measurements of secondary stable and radioactive nuclei 
in CR provide basic information necessary to probe large- 
scale Galactic properties, such as the diffusion coefficient 
and halo size, as well as mechanisms and sites of CR acceler- 
ation. Meanwhile, the accuracy of many of the nuclear cross 
sections used in CR astrophysics is far behind the accuracy of 
CR measurements of the current missions, such as Ulysses, 
ACE, and Voyager, and clearly becomes a factor restricting 
further progress. The widely used semi-phenomenological 
systematics have typical uncertainties more than ~ 50%, and 
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can sometimes be wrong by an order of magnitude (Mashnik 
2000 and references therein); this is reflected in the value of 
propagation parameters and leads to uncertainties in the in- 
terpretation. 

We have previously described a numerical model for the 
Galaxy encompassing primary and secondary CR, 7-rays and 
synchrotron radiation in a common framework (Strong et al., 
2000). Up to recently our GALPROP code handled 2 spatial 
dimensions, (R, z), together with particle momentum p. This 
was used as the basis for studies of CR reacceleration, the 
size of the halo, positrons, antiprotons, dark matter and the 
interpretation of diffuse continuum 7-rays. 

The experience gained from the original version allowed 
us to design a new version of the model, entirely rewritten 
in C++, which incorporates essential improvements over the 
older model, and in which a 3 -dimensional spatial grid can 
be employed. It is now possible to solve the full nuclear reac- 
tion network on the spatially resolved grid. We keep however 
a "2D" option since this is often a sufficient approximation 
and is much faster to compute than the full 3D case. The code 
can thus serve as a complete substitute for the conventional 
"leaky-box" or "weighted-slab" propagation models usually 
employed, with many associated advantages such as the cor- 
rect treatment of radioactive nuclei, realistic gas and source 
distributions etc. 

In this paper we show our preliminary calculations of the 
radioactive secondary ratios 26 A1/ 27 A1, 36 C1/C1, 54 Mn/Mn 
using the Los Alamos compilation of experimental cross sec- 
tions together with calculations by the code CEM2k (recog- 
nized by the nuclear physics community as among the best 
in predictive power as compared with other similar available 
codes). 

2 Model 

The GALPROP models have been described in full detail 
elsewhere (Strong and Moskalenko, 1998). The results ob- 
tained with the new version of GALPROP have been dis- 
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cussed in a recent review (Strong and Moskalenko, 2001a), 
and the most recent updates are described in Strong and Mos- 
kalenko (2001b). 

In the new version, apart from the option of a full 3D 
treatment, we have updated the cross-section code to include 
latest measurements and energy dependent fitting functions. 
The nuclear reaction network is built using the Nuclear Data 
Sheets. The beryllium and boron production was calculated 
using the authors' fits to major production cross sections 
C,N,0(p, x)Be,B. For the main channels of production of 
isotopes of Al, CI, Mn we use all experimental data avail- 
able to us from the T16 LANL compilation by Mashnik et al. 
(1998) and calculations using the improved version (Mash- 
nik and Sierk, 2000) of the Cascade-Exciton Model (Gudima 
et al., 1983) code CEM2k renormalized to the data if nec- 
essary. This code employs sophisticated microphysics via 
Monte Carlo calculations and it is difficult to use it "on-line" 
with our propagation code; other cross sections are thus cal- 
culated using the Webber et al. (1990) (wnewtr.for of 1993) 
phenomenological approximation renormalized to the data 
where it exists. For this purpose we use our internal database 
consisting of more than 2000 points collected from sources 
published in 1969-1999. (Another option is to use code 
yieldx_011000.forby Silberberg and Tsao.) For calculation 
of the total inelastic cross sections we use the latest version 
of the code CROSEC (Barashenkov and Polanski, 1994). 

The reaction network is solved starting at the heaviest nu- 
clei (i.e. 64 Ni), solving the propagation equation, computing 
all the resulting secondary source functions, and proceeding 
to the nuclei with A — 1. The procedure is repeated down to 
A = 1. In this way all secondary, tertiary etc. reactions are 
automatically accounted for. To be completely accurate for 
all isotopes, e.g. for some rare cases of 0^ -decay, the whole 
loop is repeated twice. 

3 Production cross sections 

Since the calculations with the modern nuclear codes are 
very time consuming we check the effect of the new cross 
sections only on the isotopes of Al, CI, and Mn. The ra- 
dioactive isotopes of these elements are the main astrophysi- 
cal time clocks which together with stable secondary isotopes 
allow us to probe global Galactic CR properties, in particular 
the halo size. 

For isotopes of these elements we have chosen only the 
main production channels to calculate most accurately. For 
26 Al the main progenitors are 27 Al and 28 Si, that for 27 Al is 
28 Si. For isotopes of CI the main progenitor is 56 Fe, but the 
contribution of many lighter nuclei is equally important. In 
the case of Mn, the main progenitor is 56 Fe with significant 
contribution of other isotopes of Fe, except for 55 Mn where 
only 56 Fe is important. 

The experimental data for comparison should be taken care- 
fully. The experimental technique in the past (7-spectrometry) 
did not allow for the individual partial cross section to be ex- 
tracted. Those measured represent the cross section of all the 
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Fig. 1. Production cross sections of Al isotopes. The line cod- 
ing: solid line - our adopted cross section, dashes - Webber et al. 
(1990) code (W), dash-dots - ST code. Data: T16 LANL compila- 
tion (Mashnik et al., 1998). 

nuclear reaction chains ending at the particular isotope, i.e., 
almost always cumulative yields. 

The simplest case is production of Al isotopes, however 
the abundant experimental measurements exist only for the 
natural Si -^ 26 Al reaction. Natural silicon consists of 92% of 
28 Si and the rest are isotopes 29,30 Si, 5% and 3% respectively 
(Anders and Grevesse, 1989). The measured cross section of 
the reaction " at Si(p, a;) 26 Al include also ™ at Si(p, x) 26 Si, but 
the contribution is small. 

Fig. 1 shows the experimental data for the inclusive re- 
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action nat Si(p, x) 2e Al, 27 Al(p, x) 26 Al, and 28 Si(p, x) 27 Al 
together with calculations using CEM2k and Webber et al. 
(1990) and Silberberg and Tsao (ST) codes. In case of ™ at Si, 
calculations include weighted contribution of Si isotopes. Pro- 
duction of 26 Al is calculated as the sum of 26 Al and 26 Si pro- 
duction cross sections. In this particular case we use the data 
to renormalize the 28 Si(p, x) 26 A\ cross sections calculated 
by CEM2k. In case of the reaction 28 Si(p, x) 27 Al (+ 27 Mg) 
we use a fit to the data, though the renormalized CEM2k 
model also works well above some tens of MeV. 

4 Propagation of cosmic rays 

Our preferred model for nuclei propagation is that with dif- 
fusive reacceleration. Though it has possible problems with 
secondary antiprotons and positrons (Moskalenko et al., 2001), 
it describes the spectra of nuclei and the stable secondary /pri- 
mary nuclei ratios well. We thus will use the stochastic reac- 
celeration model (SR-model) described in Moskalenko et al. 
(2001). 

As in previous work, for each halo height Zh the model 
is adjusted to fit B/C, and the source abundances at the ap- 
propriate energy adjusted to agree with the relevant observed 
stable nuclei ratios; the fluxes of the radioactive isotopes are 
then computed. In this way the uncertainty in the denomina- 
tor of the ratios is reduced. The heliospheric modulation is 
taken into account using the force-field approximation. 

Fig. 2 shows the predicted interstellar and modulated B/C 
ratio compared with observations; the reacceleration repro- 
duces the peak quite well. 

The case of 26 A1/ 27 A1 was the most uncertain giving the 
largest halo size in Strong and Moskalenko (2001a). From 
Fig. 1 it is clear that the discrepancy between the cross sec- 
tion calculations and data, which often exceeds a factor of 2, 
introduces a large error in the calculated ratio in CR. 

Figs. 3-5 show 26 A1/ 27 A1, 36 C1/C1, 54 Mn/Mn ratios cal- 
culated with the new cross sections. For this calculation we 
used the half-life times of 0.87 Myr ( 26 A1), 0.31 Myr ( 36 C1), 
0.63 Myr ( 54 Mn). The ACE data points imply a halo size of 
a few kpc. The new limits derived, 26 Al: 3.5-6 kpc, 36 C1: 
4-15 kpc, 54 Mn: 3-7 kpc, are all consistent with our limits 
derived from Be: 1 .5-6 kpc. The new limits include the error 
bars of the elemental abundance measurements from Ulysses 
(DuVernois and Thayer, 1996), to which we tune our propa- 
gated abundances. However, for isotopes of Al and Mn they 
are less important because there is only one major progenitor 
in each case. 

Fig. 6 summarizes the halo size constraints obtained in this 
analysis. These estimates are based on the four radioactive 
isotopes by requiring consistency of the calculated ratio with 
the ACE data (Yanasak et al., 2000) and taking into account 
the error bars on both prediction and data. Also shown is the 
range derived in (Strong and Moskalenko, 2001a) employing 
the Webber et al. (1990) cross section code. 

To remove uncertainty connected with elemental abundan- 
ces errors, we propose in future to derive ratios of all isotopes 
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Fig. 2. B/C ratio as calculated for a model with reacceleration. Up- 
per curve: modulated for 450 MV, lower curve: interstellar. Data: 
see Strong and Moskalenko (2001a). 



of Al, CI, Mn to the main progenitor, namely 26,27 Al/ 28 Si, 
35 - 37 Cl/ 56 Fe, and 53 ~ 55 Mn/ 56 Fe, not only the widely used 
26 A1/ 27 A1, 36 C1/C1, and 54 Mn/Mn ratios. In case of Al and 
Mn isotopes this will virtually eliminate the need to tune the 
elemental abundances. 

Some uncertainty still comes from modulation, while the 
experimental values for the ratios measured by ACE are rather 
accurate. However, because of the very flat ratio in the case 
of 54 Mn/Mn (below 1 GeV/nucleon) the modulation uncer- 
tainty is of minor importance. 

The preliminary conclusion to be drawn from all radioac- 
tive nuclei is that, at least within the context of the present 
propagation model, Zh = 4 — 6 kpc based on the ACE data 
and Ulysses elemental abundances. This is consistent with 
our previous result z h = 3 - 7 kpc (Strong and Moska- 
lenko, 2001a) and supports our previous conclusion that the 
large dispersion between the isotopes is mostly due to cross- 
section inaccuracies. 
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Fig. 3. 26 A1/ 27 A1 ratio calculated for z h = 2, 4, 6, 10 kpc (top to 
bottom). Solid curves - modulated, dashes - interstellar. Data: ACE 
- Yanasak et al. (2000), Ulysses - Simpson and Connell (1998). 
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Fig. 4. 36 C1/C1 ratio calculated for z h = 2, 4, 6, 10 kpc (top to 
bottom). Solid curves - modulated, dashes - interstellar. Data: ACE 
- Yanasak et al. (2000), Ulysses - Connell et al. (1998). 
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Fig. 5. Mn/Mn ratio calculated for z h = 2, 4, 6, 10 kpc (top to 
bottom). Solid curves -modulated, dashes - interstellar. Data: ACE 
- Yanasak et al. (2000), Ulysses - DuVernois (1997), Voyager - 
Lukasiak et al. (1997). 
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Fig. 6. Halo size limits as derived from the abundances of the four 
radioactive isotopes and ACE data. The ranges reflect errors in ratio 
measurements and source abundances. The dark shaded area indi- 
cates the range consistent with all ratios; for comparison the range 
from (Strong and Moskalenko, 2001a) is shown by light shading. 
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